Synchronization and partial synchronization of linear maps 
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We study synchronization of low-dimensional (d — 2, 3, 4) chaotic piecewise linear maps. For 
Bernoulli maps we find Lyapunov exponents and locate the synchronization transition, that numer- 
ically is found to be discontinuous (despite continuously vanishing Lyapunov exponent (s)). For tent 
maps, a limit of stability of the synchronized state is used to locate the synchronization transition 
that numerically is found to be continuous. For nonidentical tent maps at the partial synchroniza- 
tion transition, the probability distribution of the synchronization error is shown to develop highly 
singular behavior. We suggest that for nonidentical Bernoulli maps (and perhaps some other dis- 
continuous maps) partial synchronization is merely a smooth crossover rather than a well defined 
transition. More subtle analysis in the d = 4 case locates the point where the synchronized state 
becomes stable. In some cases, however, a riddled basin attractor appears, and synchronized and 
chaotic behaviors coexist. We also suggest that similar riddling of a basin of attractor might take 
place in some extended systems where it is known as stable chaos. 
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INTRODUCTION 

Recently, synchronization of chaotic dynamical sys- 
tems has been intensively studied Q . To some extent this 
is motivated by its numerous experimental realizations 
in lasers, electronic circuits or chemical reactions |2j. Of 
interest, however, are also theoretical aspects of this phe- 
nomenon. Relatively well understood is the problem of 
synchronization of identical systems. In this case it is 
known that synchronization might appear only when the 
so-called transversal, or conditional [3j Lyapunov expo- 
nents are negative [lj. Much less understood, however, 
remains the problem of partial synchronization that oc- 
curs when systems are nonidentical. In such a case 
chaotic systems do not fully synchronize which causes 
fundamental problems even with the very detection of 
synchronization. In addition to studying Lyapunov spec- 
trum, some other methods to detect this transition were 
proposed. For example it was noticed that at the partial 
synchronization of continuous dynamical systems driven 
by common noise, the probability distribution of phase 
difference changes Q. One of the difficulties in studying 
synchronization is the lack of analytical insight into this 
problem. Consequently, most of the results in this field 
are based on numerical calculations. It would be desir- 
able to study models that would be exactly solvable if 
not fully then at least with respect to some properties. 

In the present paper we study low-dimensional (d = 
2,3,4) linear maps. Such systems have the advantage 
that under some conditions their Lyapunov exponents or 
at least limits of stability of synchronized state can be 
obtained analytically. More detailed nature of synchro- 
nization transition is studied numerically. We show that 
in the case of discontinuous maps, synchronization is also 
a discontinuous transition, despite continuously vanish- 
ing (as a function of coupling strength) Lyapunov expo- 
nent (s). For continuous maps this transition is found to 
be continuous. We also study partial synchronization of 



nonidentical maps, and its relation with Lyapunov ex- 
ponents. We show that for tent maps at partial syn- 
chronization there are substantial changes in the prob- 
ability distribution of the synchronization error. But 
for Bernoulli maps partial synchronization is merely a 
smooth crossover between synchronized and nonsynchro- 
nized regimes. Finally, we examine d — 4 maps which 
might be considered as describing synchronization of two 
identical two-dimensional systems. This case requires the 
most subtle analysis, since even in the synchronized state 
the system is relatively complicated (synchronized man- 
ifold is two-dimensional in this case). Using some argu- 
ments, based on the spectral properties of Jacobian of 
this maps, we locate the point where the synchronized 
state becomes stable. In some cases, this is most likely 
a global attractor since numerical calculations show that 
synchronization transition takes place also at this point. 
However, in some cases, a basin of the synchronized at- 
tractor is riddled with basin of another chaotic attrac- 
tor ||. We suggests that such a coexistence of synchro- 
nized and nonsynchronized chaotic attractors might be a 
low-dimensional analog of the so-called stable chaos that 
was found in some spatially extended systems |(|. 



TWO-DIMENSIONAL MAPS 

Objects of our main interest are linear maps. First let 
us examine a two-dimensional map of the form 

x n +i = f((q-e)x n + (l-q)y n ) (la) 

Vn+i = f(qx„ + (1 - q)y n ) (lb) 

where f(x) is a piecewise linear function, q > 0, and 
n = 1,2... can be interpreted as a time variable. In the 
present paper we will use a generalized Bernoulli shift 
map f(x) — axmod(l) and a tent map 

ax for < x < 1/a 



fix) 



2 for 1/a < x < 1 



(2) 
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where a < 2 and in most of our numerical calculations 
we use a = 3/2. Certain aspects of dynamics of piecewise 
linear maps were already studied 0. 

First, let us consider the case of identical systems (e = 
0). Since the Bernoulli map has a constant slope (except 
a set of zero measure), Lyapunov exponents Ai and A 2 of 
the this map can be obtained from the eigenvalues of the 
Jacobian of map 

T 2 =( °* ^-d) (3) 
\a(l-q) aq ) 

Lyapunov exponents of the tent map are more difficult 
to obtain since in this case the Jacobian depends on the 
values of the arguments of the maps in fTa | -l[Tb ) . In 
particular, Jacobian takes the form only when both 
arguments are smaller than 1/a. When the argument 
of (lla|l is smaller than 1/a but the argument of Ijlbfl is 
larger, the Jacobian becomes 

tv _ ( aa «(1 -°)\ ( A \ 
L2 -{-a(l-q) -aq ) W 

Two other possibilities lead to matrices T'£ = — T 2 
(larger, smaller) and T 2 " = T 2 (both larger). To deter- 
mine Lyapunov exponents for tent map we would have 
to consider infinite product of these four matrices in or- 
der determined by the dynamics lj 1 a(l - (| 1 bll . which usually 
constitutes a formidable analytical task |8j. However, in 
our case the problem of calculating Lyapunov exponents 
simplifies greatly provided that the difference x n — y n 
is negligibly small. This is always the case in the syn- 
chronized phase and from the limits of stability of this 
phase we will be able to locate the synchronization tran- 
sition. When the difference x n — y n is small the instances 
when one of the arguments is larger and the other smaller 
than 1/a are very rare. One can thus neglect contribu- 
tions coming from matrices T 2 and T' 2 '. Consequently, 
Lyapunov exponents in this case are obtained from the 
eigenvalues of the matrix T2 only. 

Diagonalizing T2 one obtains the following Lyapunov 
exponents Ai = ln(a) with the eigenvector ei — (1, 1) 
and A2 = ln[a(2g — 1)] with e 2 = (1, —1). Let us notice, 
that e 2 is a base vector of a one dimensional transversal 
manifold on which x n 7^ y n . Moreover, A 2 = at q — 
^-1, that at the same time locates the synchronization 
transition of system (|lafl - ljlb(l both for Bernoulli and tent 
map. 

To check our results, we introduce the synchronization 
error w = (\x n — J/n|), where (. . .) denotes an average over 
time and initial conditions. Numerical results, obtained 
also for other maps described later on, were averaged over 
typically 10 3 runs with different initial conditions. Each 
run has a length 10 7 steps with 10 6 steps discarded as 
initial transient. In some cases, where fluctuations were 
strong or relaxation time was very long we made runs of 
10 8 steps with 10 7 steps discarded. 



Results of the numerical calculations for a — 3/2 are 
shown in Fig. ^ They confirm that our system i|la|) - 
(flb)) synchronizes at q = = 5/6(= 0.833 . . .) both 
for Bernoulli and tent maps. Let us notice, however, a 
difference in the behavior of these two maps. Namely, 
while for tent map w vanishes continuously to zero at 
q = 5/6 it has a discontinuous jump for Bernoulli map, 
although A2 vanishes continuously to zero in this case. 

Now, let us consider the case of nonidentical systems, 
i.e., e > 0. In this case the Jacobian for the Bernoulli 
map has the form 

T 2 = ( °(« - e \ a{l ^ ) (5) 
\a(l-q) aq J 

Its diagonalization gives the following Lyapunov expo- 
nents 

Ai, 2 = ln[| (2q-e± y/? + 4(1 - a) 2 )] (6) 

For a = 3/2 elementary calculations show that now 
A 2 < for 



As before Ai remains positive. 

Since our systems are now nonidentical, the difference 
x n —yn never vanishes and neglecting matrices T' 2 and T 2 
is this time only an approximation. But we expect that 
for small e and in the partially synchronized state this 
approximation should be relatively accurate. Thus, we 
assume that Eq. © gives approximate Lyapunov expo- 
nents at least inside partially synchronized phase. Hence 
Eq. J7J locates the partial synchronization transition also 
for tent map. 

Numerical calculations of w as a function of q are 
shown in Fig. ^ As expected, w remains positive even 
for q < q c although for tent map there is a clear change 
of behavior around q c . For Bernoulli map w seems to be 
a smooth function of q even at q c . 

Recently, it was shown that partial synchronization in 
some continuous dynamical systems manifests through 
changes in the probability distribution of the phase differ- 
ence y . To check whether this is more general property 
we calculated probability distribution of the synchroniza- 
tion error P(\x n — y n \) and the results are shown in Fig-El 
One can see that in the nonsynchronized phase of both 
tent and Bernoulli map P(\x n — y n \) has a rather broad 
and smooth distribution. However, for tent map at q = q c 
this distribution becomes very complex. This complex- 
ity persists also in the synchronized phase (q < q c ). On 
the other hand, there is no noticeable difference at q = q c 
for the Bernoulli map, although it develops a pronounced 
peak in the synchronized phase. Taking into account also 
the smooth behavior of w at q = q c it is legitimate to ask 
whether partial synchronization in the case of Bernoulli 
map is a sharp and well defined transition or it is rather 
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FIG. 1: The synchronization error toasa function of q for the 
two-dimensional map l|la^ - ljlb(l with a — 3/2. Inset shows the 
behavior of w for tent map in the vicinity of the transitions 
for e = (x) and e = 0.01 (□). For e = 0.01 A 2 from Eq. © 
vanishes at q = q c — 0.83587 . . . that is indicated by a small 
arrow. 
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FIG. 2: The nonnormalized probability distribution P(\x„ — 
y n \) for tent (a) and Bernoulli (b) map (a = 3/2 and q c = 
0.83587...). 

a smooth crossover. In this context it would be desirable 
to examine a map that continuously interpolates between 
these two maps as for example a skewed tent map. 

THREE-DIMENSIONAL MAPS 

As a next example let us consider a three-dimensional 
map of the form 

x n +i = f((q -e)x n + 0.5(1 -q)(y n +z„)) (8a) 
Un+i = f(qy n + 0.5(1 - q)(y n + x n )) (8b) 
z n +i = f(qz n + 0.5(1 -q)(x n + y n )) (8c) 



As before, Lyapunov exponents for the Bernoulli map are 
obtained from the eigenvalues of the following matrix 

/ a(g-e) a(l - q)/2 a(l - q)/2 \ 
T 3 = a(l-g)/2 aq a(l - g)/2 (9) 
\a(l-g)/2 a(l-g)/2 aq J 

And similarly to the d = 2 case, for e = and in the 
synchronized phase, matrix T3 contains also Lyapunov 
exponents of tent map (in the synchronized state all three 
variables are the same). 

For e = we obtain that the largest exponent Ai 
is always positive. Two other exponents are degener- 
ate A2 = A3, their eigenvectors are the base of a two- 
dimensional transversal manifold, and for a — 3/2 they 
vanish at q = 7/9(= 0.777...). Numerical calculation 
of the synchronization error show (see Fig. that at 
this point, as expected, synchronization transition takes 
place. Similarly to d = 2 map, for Bernoulli map this 
transition is discontinuous. 

For nonidentical systems (e = 0.1) we find that A3 
vanishes at q = q c = gTjjf^ = 0.82717 . . .. On similar 
grounds as in the d = 2 case, we consider it also as an 
approximate location of the partial synchronization tran- 
sition in the tent map. Figure confirms this argument 
since one can see that synchronization error uu for tent 
map indeed changes its behavior very close to this point. 
For Bernoulli map w seems to be again a smooth function 
around q c . But there is yet another transversal exponent, 
namely A2 (with the eigenvector e2=(0,l,-l)) that van- 
ishes as in the e = case i.e., at q = 7/9. At this point 
complete synchronization of y and z systems takes place 
(they are identical). At this point w has a jump but only 
for Bernoulli map. 

We also examined the probability distribution P(x n — 
y n ). Our results (not shown) indicate that for tent map 
P(\x n — Vn\) nas a complex structure in the synchronized 
phase (q < q c ). On the other hand no appreciable change 
is seen for Bernoulli map which is similar to the d = 
2 case that we already discussed. It suggests that for 
discontinuous maps there are almost no indications of a 
qualitative change at partial synchronization (except of 
the vanishing Lyapunov exponent) and most likely it is 
only a smooth crossover. 



FOUR-DIMENSIONAL MAPS 

As a final example let us consider a four-dimensional 
map of the form 

x n +i = f{0.9x n + qy n + 0.lz n ) (10a) 

2M+i = /(0.9y„ + qx n + 0.1u„) (10b) 

z n +l = /(0.9z„ + qu n + 0.lx n ) (10c) 

u n+ i = f(0.9u n + qz n + 0.ly n ) (lOd) 
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FIG. 3: The synchronization error w as a function of q for 
the three-dimensional map l|8a fl -l|8c [l . For e = 0.1 A3 vanishes 
at q — q c = 0.82717 . . . indicated by an arrow (a = 3/2). 

This map describes a coupling of four identical systems 
x, y, z, and u. Increasing q increases the coupling x — y 
and z — u and for sufficiently large q the system enters a 
synchronized state were x n = y n and z n = u n . Numerical 
coefficients in this map are chosen in such a way that 
subsystems x — z and y — u remain nonsynchronized and 
chaotic. 

Lyapunov exponents for the Bernoulli map are ob- 
tained from the eigenvalues of the following matrix 

/ 0.9a aq 0.1a \ 

aq 0.9a 0.1a 

4 0.1a 0.9a aq K ' 

\ 0.1a aq 0.9a / 

Diagonalizing T4 we find that two Lyapunov exponents 
Ai = log[a(l + q)} and A 2 = log[a(0.8 + g)] are always 
positive. Corresponding eigenvectors ei = (1, 1, 1, 1) and 
62 = (1, 1, — 1, — 1) are the base of the synchronized sub- 
space (on which x n — y n and z n — u n ). There are 
also two transversal exponents A3 = log[a(l — q)] and 
A4 = log[a(0.8 — q)] with eigenvectors e 3 = (—1, 1, —1, 1) 
and e4 = (1,-1,-1,1). For a = 3/2 one can see that 
A 3 = at q = 1/3 and A 4 = at q = 2/15(= 0.1333 . . .). 
One can thus expect that the system synchronizes when 
the largest transversal exponent vanishes i.e., at q = 1/3. 
Indeed, such a behavior is seen for the Bernoulli map 
(Fig. @J. In this case one can also see a jump in w at 
q = 2/15 that is a trace of vanishing of A4. However, for 
tent map a different behavior is seen, and synchroniza- 
tion takes place around q = 0.23. 

To explain such a behavior we have to notice that the 
synchronized manifold is in this case two dimensional and 
in the synchronized state the system is described by two 
variables e.g., x n (= y n ) and z n (— u n ). It means that in 
the synchronized state the difference x n — z n does not 
vanish. Consequently, arguments of maps (|10afl - l|10d|) 
might differ and we have to take into account that in 



some cases Jacobian of the map has a different form than 
T4. For example, when arguments of the first two maps 
p()a|l - l|10b|l are both smaller than 1/a and arguments 
of maps are both larger than 1/a, then the 

Jacobian is given by the following matrix 

/ 0.9a aq 0.1a \ 

, = aq 0.9a 0.1a , . 

4 -0.1a -0.9a -aq K ' 

\ -0.1a -aq -0.9a / 

In the following, similarly to the d = 2 and 3 cases, we 
will analyse the stability of the synchronized state. When 
differences x n — y n and z n — u n are negligibly small then 
arguments of maps l|10a|) and IjlObfl are almost the same, 
and arguments of maps (|10c|) and ijlOdfl are also almost 
the same. Thus, the situation is analogous to the d = 2 
example and there are only four possible forms of the 
Jacobian given by the matrices T4, T 4 , T 4 = — T 4 , 
and T'l' = -T 4 . However, since T 4 and T^' do not 
commute with T 4 and T 4 finding Lyapunov exponents 
of our map even under such a simplifying assumption is 
very difficult. But a closer look at the spectral properties 
of these matrices reveals certain interesting properties 
(with respect to their spectrum, it is of course sufficient 
to examine for example only matrices T4 and T 4 ). Cal- 
culating eigenvalues of T 4 one finds the following Lya- 
punov exponents: \[ = \' 2 = \xi[a\J (1 + q)(0.8 + q)], 
A 3 = A4 = ln[ay (1 — q)(0.8 — q)\. These Lyapunov ex- 
ponents are associated only with matrix T 4 and they are 
not Lyapunov exponents of map (|10afl - l|10dl) . Exponents 
Xi and X' 2 are positive and their eigenvectors are within 
the synchronized manifold. Eigenvectors of exponents 
A 3 and A 4 are within the transversal manifold. These 
exponents change sign at q = q c = 0.9 — \/0.01 + a -2 , 
that for a = 3/2 equals 0.225875 . . ., i.e., the value that 
is very close to the actual transition in the tent map as 
obtained numerically (Fig0J. However, for such a value 
of q only one of the transversal exponents of T 4 is neg- 
ative (A 4 ) and the other (A3) remains positive. But for 
q = 0.9 — VO.01 + a~ 2 and interesting property of matrix 
T 4 holds namely A3 + A4 = 0. It suggests that indeed, at 
this value a synchronization transition of the tent map 
might take place, however a closer look into the dynamics 
of our map HlOafl - ljlOdfl ^ s needed. 

Let us analyse the dynamics of our system for q > 

0. 9 — VO.01 + a~ 2 as projected on the transversal man- 
ifold. Each state vector y n , z„, u n ) is decomposed 
in the base of eigenvectors of T4 as (x n , y n , z n , u n ) = 
ciei + C2e2 + 0363 + 0464, and then we look at the evolu- 
tion of transversal coordinates C3 and C4. A typical tra- 
jectory is shown in Fig. [S] One can see that after an ini- 
tial transient, the system locates essentially on the axes 
of eigenvectors and gradually approaches the point (0,0) 

1. e., a synchronized state, swapping from time to time its 
coordinates. Actually, such a behavior can be deduced 
from Eqs. I|10a|) - (|10d|) . First, we express coordinates in 
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FIG. 4: The synchronization error w as a function of q for 
the four-dimensional map l|10a|l - l|10d|l . We also measured the 
difference (\z„ — u n \) and found that within statistical errors, 
as expected, it equals to w. An arrow indicates a value q = 
0.225875 . . .. 
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FIG. 5: A typical trajectory in a transversal subspace 
( | C3 1 , | C4 1 ) for q — 0.25 and a = 3/2 (only absolute values 
of coordinates C3,C4 are plotted). Arrows indicate the action 
of Lyapunov exponents of T4. Due to swapping of coordi- 
nates (see text) the system evolves toward the synchronized 
state (0,0). 



the transversal subspace as C3 — (—x n +y n — z n +u n )/4 = 
(-r5i-5 2 )/4 andc 4 = (x n -y n - z n + u n )/A = (5i — <5 2 )/4, 
where 5% = x n — y n and ^2 = z n — u n . Since 61,62 
are small it implies that the system evolves 'in pairs'. 
Namely, either both arguments of maps (|10a|) - ljl0b|l are 
smaller or both are larger than 1/a and either the first 
or the second definition of the tent map J2J is applied to 
this pair. The same applies for the second pair of maps 
(I10c|) - (ll0d|) . The swap of coordinates appears when two 
pairs evolve according to different rules. Indeed, let us 
assume for example that the first pair of arguments is 
smaller and the second pair is larger than 1/a. Using 
the evolution equations (|10a|) - (|10d[) we find that initial 
coordinates ((—Si — <5 2 )/4, (<5i — <52)/4) are transformed 
into (-a(0.8 - q)(5i - 5 2 )/4 } a(l - q){S 1 + <5 2 )/4). Up 
to the sign and the multiplicative factor of the order of 
unity this is indeed the swap of coordinates. Since for 
q > 0.9 — VO.01 + a~ 2 the sum of transversal exponents 
A3+A4 < 0, shrinking along C4 axis prevails over repelling 
along C3 axis. Without swapping the system would drift 
to infinity along C3 axis. But swapping turns the sys- 
tem back almost on the C4 axis where it is subjected 
(mainly) to shrinking. Of course, swapping can also put 
the system again on the C3 axis where it is repelled, but 
sooner or later it will return on the C4 axis, where it 
will be again subjected to shrinking. Swapping of co- 
ordinates is actually described by matrices T4 or T4. 
For q > 0.9 — VO.01 + a~ 2 these matrices have negative 
transversal exponents and their action in the transver- 
sal manifold can only shrink the system toward the ori- 
gin (0,0). As a result of this combined action (swapping 
and shrinking faster than repelling) , the system must end 
up in the synchronized state. On the other hand for 
q < 0.9 — \/0.01 + a~ 2 both mechanisms which attracted 



the system toward the synchronized state fail. Repelling 
along C3 axis prevails over shrinking along along C4 axis 
and transversal exponents of T4 and T4 are positive. 
As a result the system remains in the nonsynchronizcd 
state. The value q = 0.9 - V0.01 + a- 2 = 0.225875 . . . 
(for a = 3/2) is indicated in Fig. ^ with an arrow and is 
in a very good agreement with numerical calculations. 

There is, however, a silent assumption in the above 
considerations, namely that on the way to the synchro- 
nized state the system will not fall into the basin of an- 
other attractor that might possibly exist. Actually, it 
is already known that basins of synchronized chaotic at- 
tractors are sometimes riddled with basins of some other 
attractors, even when transversal Lyapunov exponents 
are negative [f|. As a result, only some initial condi- 
tions will lead the system to the synchronized state, and 
the other will lead to the other nonsynchronized attrac- 
tor. And indeed, we found that for a = 1.35 such a 
scenario most likely takes place (Fig. |SJ). For a = 1.35 
we obtain that the synchronized state becomes stable for 
q > 0.9 - V0.01 + a- 2 = 0.1525397. . .. Although the 
synchronization error w does not vanish at this point, we 
observed that for q > 0.154 there is a finite fraction P 
of initial conditions that leads relatively fast to the syn- 
chronized state. Remaining initial conditions lead the 
system to another chaotic attractor. As a result, w that 
is an average over runs with different initial conditions, 
remains positive. 

To show that a — 1.35 and a = 1.5 cases are qualita- 
tively different we plot attractors at their critical points 
q = q c = 0.9 — VO.01 + a' 2 as projected on the transver- 
sal manifold (03,04) (Fig. 0). These plots are obtained 
from several independent runs after discarding 10 8 steps 
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FIG. 7: Critical point attractors as projected on the transver- 
sal manifold (03,04). 
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FIG. 6: The synchronization error w and the probability P of 
reaching the synchronized state as a function of q for the four- 
dimensional map (llUal - dlOdll with a — 1.35. Arrow indicates 
a value q — 0.1525397 . . ., where a synchronized state becomes 
stable. 



and recording next 3 • 10 4 steps. We noticed, that when 
q is smaller than q c no appreciable change is seen in the 
structure of these attractors. However, when q exceeds 
even slightly the critical value important differences takes 
place. Namely for a = 1.5 the attractor shrinks to the 
origin (0,0). For a = 1.35 the origin also becomes an 
attractor but the nonsynchronized attractors persist as 
well, with only small changes. In Fig. probability P of 
reaching a synchronized state seems to be positive only 
for q > 0.154, which might be a consequence of the nu- 
merical method. Namely we considered that a system 
reached a synchronized state if after 10 6 steps its coor- 
dinates satisfy \x n — y n \ + \z n — u n \ < 10~ 5 . But close 
to the critical poind q c an approach to the synchronized 
state is very slow and we expect that when larger num- 
ber of iteration steps were used, we would find positive 
P for any q > q c . Let us also notice that the fact that 
for a = 1.5 at the critical point the attractor is different 
than the origin implies a finite value of synchronization 
error w at this point i.e., the synchronization transition 
is most likely discontinuous in this case. 

Finally, let us notice that in synchronization of some 
spatially extended systems, like Coupled-Map Lattices, 
it was observed that in presence of negative transver- 
sal Lyapunov exponents nonsynchronized state might ex- 
ist A similar behavior appears in some other ex- 
tended systems and is known as stable chaos 0. In our 
opinion, the behavior of our system for a — 1.35 bears 
some similarity to this effect. Although we did not cal- 
culate the transversal Lyapunov exponents explicitly, it 
seems very likely that for q > q c they are negative (since 
the synchronized state is stable, and that is also con- 
firmed numerically because a finite fraction of initial con- 
ditions leads to this state). Nevertheless, a finite fraction 



of initial conditions lead the system to a nonsynchronized 
state. An important difference is that in extended sys- 
tems exhibiting stable chaos a nonsynchronized state is a 
generic state of the dynamics, that is reached from almost 
any initial condition. Another property of such extended 
systems is that a nonsynchronized state after a certain 
threshold looses its stability but no such a threshold was 
found in our 4-dimensional map. But it might be that 
these features are consequences of high dimensionality of 
dynamics of an extended system. 



CONCLUSIONS 

In conclusion, we studied synchronization of low- 
dimensional linear chaotic maps. Even for such simple 
maps finding Lyapunov exponents might be very diffi- 
cult. But as we have shown, the problem significantly 
simplifies in the synchronized phase which enabled us to 
locate exactly the synchronization transitions. Our re- 
sults show qualitative differences between the behavior of 
Bernoulli and tent map concerning the type of the syn- 
chronization transition (discontinuous for Bernoulli map, 
and continuous for tent map) , or even the location of the 
transition (in the d = 4 case). Morever, very different 
behavior is seen when maps are nonidentical and only 
partial synchronization takes place. It might be that it 
is a discontinuity of the Bernoulli map that is responsi- 
ble for some of these differences. Noting the similarity 
to the behavior of some extended systems, it would be 
interesting to examine whether the origin of the stable 
chaos indeed can be relateded with riddled basins of at- 
tractors in low-dimensional systems like those examined 
in the present work. 
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